跳到主要内容

Gram-Schmidt 正交化

阐述​

  • q1q_1 由 a1a_1 归一化得到
  • q2q_2 由 a2−q1q1∗a2a_2-q_1q_1^*a_2 归一化得到
  • q3q_3 由 a3−q1q1∗a3−q2q2∗a3a_3-q_1q_1^*a_3-q_2q_2^*a_3 归一化得到
  • …… 

从这些就可以算出 RR 的系数。

算法:

这大约需要 2mn22mn^2 次浮点运算。但是,问题在于这样相当于不断用 Rj−1R_j^{-1} 来右乘 AA,其中有可能 rjj≈0r_{jj}\approx 0,这样会导致 RR 是病态的。

所以,可以从另一个方面着手,用正交矩阵来左乘 AA,这样可以保证 κ(Qj)=1\kappa(Q_j)=1.

实例​

性质​

可以用于 QR 分解。

相关内容​

参考文献​